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Abstract 

In this article we investigate averaging properties of fully nonlinear PDEs in bounded domains 
with oscillatory Neumann boundary data. The oscillation is periodic and is present both in the 
operator and in the Neumann data. Our main result states that, when the domain does not have flat 
boundary parts and when the homogenized operator is rotation invariant, the solutions uniformly 
converge to the homogenized solution solving a Neumann boundary problem. Furthermore we 
show that the homogenized Neumann data is continuous with respect to the normal direction of 
the boundary. Our result is the nonlinear version of the classical result in [4] for divergence-form 
operators with co-normal boundary data. The main ingredients in our analysis are the estimate 
on the oscillation on the solutions in half-spaces (Theorem |3.1[ ), and the estimate on the mode 
of convergence of the solutions as the normal of the half-space varies over irrational directions 



(Theorem 4.1 1 



1 Introduction 



Let us consider a bounded domain in R" containing the closed unit ball K = {x : \x\ < 1} (see 
Figure 1). Let g : M" — [1,2] be a Holder continuous function which is periodic with respect to the 
orthonomal basis {(ei, e„)} of M". More precisely, g satisfies 

g{x + ei) = g{x) for i = 1, n and g G C'^(M") for some < /? < 1. 

With r2, K and g as given above, we are interested in the limiting behavior of the following problem: 

' F{D^u'', f ) = inn-K, 



3(f) 



on (90. 



Here v = is the outward normal vector at a; S 9f2 and F is a uniformly elliptic, fully nonlinear 
operator (see section 2 for the precise definitions and conditions on F). The set K is introduced to 
avoid discussions of compatibility conditions on the Neumann boundary data. The operators under 
discussion include, for example, the divergence or non-divergence form operator 



'^i,j^ij{~)dxiXj'U — 0, 



(1) 
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Figure 1 



where aij(x) are Holder continuous and XInxn < {0'ij)ix) ^ ^Inxn for positive constants A and A. 
Let us define 

5""^ = {i^ e M" : = 1} and M" := {M : n x n symmetric matrix }. 

Due to [5], there exists a uniformly elliptic operator F : Ai" M such that any solution of 
F{D^u', -) = in with fixed boundary data on Qfi converges uniformly to the solution of F{D'^u) — 



in n (see Theorem 2.12 1. Therefore the question under investigation is whether the oscillatory boundary 
data changes the averaging behavior of w^. As we will discuss below, several difficulties arise in answering 
this question due to the nonlinear (or non-divergence form) nature of the problem as well as the geometry 
of the domain. 

Let us state our main result. First we introduce the following "cell problem" for given v € S"~^ 
and A e M. 



' FiD'^u", f ) = in{-l<x-v-\< 0}; 

~ 1 on r^i :— {x ■ V = \ — 1\; 

d^u^^gi^) ouTn := {x ■ ly ^ X}. 

[ and, {P)£ and {P)e,u.\.a be as given above. 



Theorem 1.1. Let F : A^" - 

(a) [Theorem 3.1] For given normal direction v G S"^^^ — MZ" and any A G M there exists a unique 
homogenized neumann data ^{v) for solutions of {Pe,v.\). 

(b) [Theorem 5.1] Suppose that F{M) only depends on the eigenvalues of M G A^". Then there 
exists a continuous function fL{v) : 5""^ — > M., given as the continuous extension of fi{v) over 
V G S"^^^ — MZ", such that the following holds: 

Suppose n is a bounded domain in K" such that dfl is and does not contain any flat part. 
Let u"^ solve {P)e- Then converges locally uniformly to u, which is the unique solution of the 
homogenized problem 



(P) 



F{D^u) =0 infl-K; 
u = 1 on K] 

duU — fi{i') on dfl. 
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Remarks 1.2. 1. The assumption on F in Theorem \ (b) seems to he necessary condition to 
achieve the continuity of the homogenized slope p-iiy), at least from numerical experiments (see 
JJBl)- Note that this assumption is equivalent to the rotation and reflection invariance of F. For 
more discussion on the role of the assumption in the analysis, see the discussion in section 4- 

2. We point out that the restriction on the geometry of dQ is necessary: in fact, if dfl is locally 
T = {x ■ v = X} with V G Z" and \ ^ Q, then the (normalized) distribution of boundary data ) 
on r changes a lot as e — > 0. Consequently in this case there is no unique limit of . 

oDiscussions on previous results 

Our problem is classical for the case of uniformly elliptic, divergence-form equations 

- V • (A(-)Vun = in 17, (2) 
e 

with the co-normal boundary data 

ly- {A{-)Vu')ix) = g{-) on dn. (3) 

£ e 

For ([2|-([3]), a corresponding result to Theorem 1 1 . 1 1 was proved in the classical paper of Bensoussan, 
Lions and Papanicolau [1] with explicit integral formula for the limiting operator as well as the limiting 
boundary data. Recently, a corresponding result was shown for systems of divergence type operators 
with oscillatory Dirichlet boundary data in convex domains, by Gerard- Varet and Masmoudi [14] (also 
see [IT]). 

For nonlinear or non-divergence type operators, or even for linear operators with oscillatory Neu- 
mann boundary data that is not co-normal, most available results concern half-space type domains 
whose boundary goes through the origin. In [19 , Tanaka considered some model problems in half-space 
whose boundary is parallel to the axes of the periodicity by purely probabilistic methods. In T] Arisawa 
studied special cases of problems in oscillatory domains near half spaces going through the origin, using 
viscosity solutions as well as stochastic control theory. Generalizing the results of [T], Barles, Da Lio and 
Souganidis 3J studied the problem for operators with oscillating coefficients, in half-space type domains 
whose boundary is parallel to the axes of periodicity, with a series of assumptions which guarantee the 
existence of approximate corrector. In |10| the continuity property of the averaged Neumann boundary 
data in half-spaces, with respect to the normal direction, was studied in the case of homogeneous oper- 
ator F. Recently, Barles and Mironescu [5j showed a corresponding result to (3j for oscillatory Dirichlet 
boundary data. 

o Main ingredients and challenges 

The main steps in extending aforementioned results from half-space type domains to general domains 
are the following. 

The plan is to use these solutions in strip domains to approximate those in general domains. For the 
stability of such approximation one requires the distribution of g to be invariant on Fjv regardless of 
the choice of A, which is the case if v is not a multiple of vectors in Z": we call such vectors irrational. 
For irrational we show that there is a unique linear profile the solutions converges to, by proving an 
oscillation lemma as well as using a quantitative version of Weyl's equi-distribution theorem. 

Secondly, to establish sufficient stability to address the general domains, some estimate on the 
convergence mode of e- solutions in the strip domain, in terms of the variation of i', is necessary. 
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Establishing this estimate is our second, and most chahenging, main step (Theorem 4.1). For problems 
with divergence-form structure, such estimates were obtained in [3] and [H] (also see [T7]) by means 
of integral formulas. In our setting we must proceed by maximum principle-type arguments, which 
requires careful perturbation of the boundary data as well as construction of delicate barriers which 
describes the averaging behavior of solutions up to the Neumann boundary. The proof of Theorem 4.1 
is based on the observation that, in "meso-scopic" scale, the distribution of g on two hyperplanes with 
close irrational normals are similar (see section 4.1. for a heuristic discussion of this fact). We adopt 
a multi-scale homogenization argument to address separately the effects on the solution caused by (a) 
microscopic oscillation of g near the Neumann boundary (b) the difference in normal directions and (c) 
the oscillation present in the operator F . 

We mention that a significant difficulty arises due to the presence of the ^-dependence in F in 
addition to the oscillations in g. To get around this difficulty, we use localization arguments as well 
as the existing homogenization results (e.g. [9] and 8J) to show homogenization occurs away from the 
Neumann boundary. Such strategy works since near the Neumann boundary the oscillation of the first 
derivative {g) dominates the behavior of solutions. 

o Outline of the paper 

In section 2 we introduce some notations as well as preliminary results which will be used in the 
rest of the paper. In section 3 we first study the averaging properties of the operator in the strip 
(half-space type) domains, to show that if the hypersurface is normal to an irrational direction then 
there is a unique homogenized Neumann data in the limit e — > 0. Besides the complications that the 
inhomogeities in the operator F cause, the proof of averaging phenomena in this setting is due to the 
Weyl's distribution theorem, whose quantitative version that we need is borrowed from [TD] and is stated 



in Theorem 2.2 In section 4 we prove the main estimate on the mode of convergence of e-solutions in 
the strip domain as the normal v varies around a reference direction. After a heuristic description of 
the proof, we prove a series of lemmas which describe the behavior of solutions in different regions of 
the strip domain, divided in terms of the distance to the Neumann boundary. Based on the lemmas 
then we are able to construct suitable barriers to obtain the desired estimate (Theorem 4.1). Lastly in 
section 5 we prove Theorem 1 1 . 1 1 using the aforementioned estimate, the non-flat geometry of dVl, and 
the stability of viscosity solutions. 

2 Preliminaries 

2.1 Rational and Irrational directions 

We first introduce the following categories of normal directions, following [lOj . 
Definition 2.1. (i) v E S''^^^ is a rational direction if v E MZ" . 



(ii) V E S" ^ is an irrational direction if v is not a rational direction. 

Hyperplanes with irrational directions represent surfaces on which the boundary data does not 
change too much when one translates it in the normal direction. The following lemma quantifies this 
fact. The proof is a straightforward application of Weyl's equi-distribution theorem. 

Lemma 2.2 (Lemma 2.7, ^Q\). For v e M" and xq G M", let us define H{xo) ■= {x : {x - xq) ■ = 0}. 
Then the following is true for < e < 1 ; 
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(i) Suppose that v is a rational direction. Then there exists a constant > depending only on v 
such that, for any x e H{xq), there is y G H{xq) satisfying 

\x-y\<M^; y-XQ&T". 

(ii) Suppose that v is an irrational direction. Then there exists a mode of continuity Wi, : [0, 1) — )■ M+ 
and a dimensional constant M > such that the following is true: for any x G H{xq), there exists 
yeW such that 

|a; - 2/1 < M£-9/^°; y - S Z" 

and 

d.ist{y,H{xQ)) <uj^{e). (4) 
(Hi) If u is an irrational direction, then for any a; e M" and 5 >Q, there exists y G H{xo) such that 

\x-y\<5 mod Z". 

Remarks 2.3. The power —9/10 appearing in (ii) of above lemma is arbitrarily chosen, and can be 
replaced with any < a < 1. The mode of continuity cj^, of course, depends on the choice of a. 

2.2 Well-posedness and regulEirity results 

Next we introduce the properties of the operator F. As before, let A4" denote the normed space of 
symmetric nx n matrices. In this paper we assume that the function F{M, y) : M.^ x R" M satisfies 
the following conditions: 

(Fl) [Homogeneity] F{tM,x) = tF{M,x) for any M € A^", x G M" and t > 0. In particular it follows 
that F(0, x) = 0. 

(F2) [Lipschitz continuity] F is locally Lipschitz continuous in M.^ x R", and there exists a constant 
C > such that for all x,y and M,N eM"" 

\F{M,x) - F{M,y)\ < C{\x - y\{l + \\M\\ + \\N\\) + \\M - N\\). 

(F3) [Uniform EUipticity] There exists constants < A < A such that 

Ajj II < F{M, x) - F{M + N,x)< A\\N\\ for any M,N ^M^ and N >Q. (5) 

(F4) [Periodicity] For any M e and x e M", 

F{M, x) = F{M, x + z)ifzG Z". 

A typical example of F satisfying (F1)-(F4) is the non-divergence type operator 

F{D^u,y) = -T,ijaij{y)dxix,u, (6) 

where aij : M" ^> M is periodic, Lipschitz continuous and A = [aij] satisfies X{Id)nxn l£ A < A(/rf)„xn- 
Another example is the Bellman-Issacs operators arising from stochastic optimal control and differential 
games 

F{D^u,y) = inf sup{£"'^u}, 
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where C^^ is a two-parameter family of uniformly elliptic operators of the form ([6]). We refer to [TT] 
for further examples of F. 

Let 51 and K be as given before, and let /(y, v) : K" x 5""^ — >■ M be a continuous function. To 
incorporate both {P)e and the homogenized problem (P), let us introduce a definition of viscosity 
solutions for the following problem: 



' F{D^u,y) = in fl-K; 
V ■ Du — f{x, v) on dfl; 
u = 1 on K, 



where v = is the outward normal at the boundary point x G dfl. See [12] for a game-interpretation 
of (-P)/. The following definition is equivalent to the ones given in [TT] : 

Definition 2.4. (a) An upper semi- continuous function w : f2 — > M is a viscosity subsolution of (P) f 

if u cannot cross from below any (j) G C^(il) H C^{Q.) which satisfies 

F{D'^(t),x) > m n- K: v ■ Dcj) > f{x,v) on d^; (j> > 1 on K. 

(b) A lower semi- continuous function u : fl ^ M is a viscosity supersolution of (P) f if u cannot cross 
from above any Lp G C'^{VL) C\ C^{il) which satisfies 

F{D'^ip, x) < in - K; v ■ Dtp < f{x, v) on (9f2; ip <l on K. 

(c) u is a viscosity solution of {P)f if u is both a viscosity sub- and supersolution of {P)f. 

The existence and uniqueness of viscosity solutions in bounded domains are consequences of the 
following comparison principle. 

Tiieorem 2.5. [Section V, 1 161] Let T, be a domain in M" whose boundary consists ofT^ and Td, and 
let F satisfy (Fl) — {F4). Suppose u, u : E — > M are bounded and continuous. If u and v satisfy 

(a) F{D'^u,x) < < F{D'^v,x) in VL; 



(b) u < V on r 



D, 



(C) t<fi^)<to^^N. 

where f{x) : M" — > M is continuous, and (a) and (c) should be interpreted in the viscosity sense. Then 
u < V in n. 

Next we state regularity results for solutions of F{D^u, x) — which will be used later. Let us recall 
the Pucci extremal operators: for a symmetric n x n matrix M, we define 

V+{M) := A(I],.>oe,) + A(Ee.<oe.) 

and 

V-{M) := A(Ee,>oe,) + A(Ee,<oeO, 

where {e^j's are the eigenvalues of M. The operator —V^ then satisfies the assumptions (F1)-(F4). 
Moreover note that, due to ([5]), we have 

- -p+(Af - N) < F{M, x) - F{N, x) < -V^ (M - N) for any M, N e M''\ (7) 

In particular, the difference w of two viscosity solutions oi F{D^u, a;) = satisfies both —V^{D'^w) < 
and —V~{D^w) > in the viscosity sense. 
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Theorem 2.6 (Proposition 4.10, [5]: modified for our setting). Let u be a viscosity solution of 

-V+{D'^u)<Q, -V-{D'^u)>0 (8) 
in a domain Q.. Then for any < a < 1 and a compact subset Q! ofil, we have 

liu||c°(n') < C'd""||M||Loo(o) < oo 
with d = d(ri' , dn) and C > depending only on n, A, A. 

Theorem 2.7 (Theorem 8.3, [6]). Let F satisfy (Fl) — (FA), and u be a continuous function in Bi 
solving 

F{D^u,x) = m Bi(0). 
Then there exists < a < 1 and C > 0, depending only on A, A and n, such that 

ll"llci.°(i3i/2(0)) < C'II'«IIl~(Bi)- 

Theorem 2.8 (Theorem 8.2, [IB]: modified for our setting). Let 

B+ := {\x\ < 1} n {a; • e„ > 0} and T = {|a;| < 1} n {a; • e„ = 0}. 

For given bounded function g : — > M, the following is true for u satisfying 

-V+{D^u)<0, -V-{D^u)>Q m B+; 
V ■ Du = g onT. 



(a) u £ C"{B^^^) for < a = a{72,X,A) < 1 given in Theorem 



2.7 with the estimate 



(b) Lf, in addition, g £ C^{B^) for some < (3 < 1 and u satisfies F{D^u,x) — 0, then u € 
Ci'"(B+2), where a = min(aoj/3) with ao given in (a). Moreover, we have the estimate 

where C is a constant depending only on n, A, A and a. 

2.3 Localization 

Here we discuss how to locahze solutions of Neumann-boundary problems. For given i' E 5"^^, let us 
define 

n^{0):={x:~l<x-iy<0}, Tjv = {a; • = 0}, Ti, = {a; • = -1}. 
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Lemma 2.9. For given R,e > 0, suppose u : Tli,{0) — > M solves the following: 

' -V+{D^u)<0, -V-{D^u)>0 in n^(0) n {|a;| < i?}, 
d,yU — g{x) on Fat n {|a;| < R\; 

u^f{x) onTDn{\x\ < R}, 

\u\{x)<R'^-^ on\x\=R. 
If 1/(2;) I, \g{x)\ < 5, then we have 

\u{x)\ < 25+ 16(n- in n^(0) n {\x\ < 1}. 

A 

Proof. Without loss of generality, let us set u — e^, x' — (xi, ...,0:^-1), .and consider the function 

^{x) = {S + 4(n - l)^i?-^)(x„ + 2) + ^{\x'\' ^{n- l)^((x„ + 1)^ - 1)). 

Then -V+iD'^ip) = 0, and > (5 on {a;„ = -1}, ip > R^-^ on = R, and d^^ip > 5 on {x,, = 0}. 
Therefore by comparison principle (Theorem 2.5) we have 

u<(pin n {|x| < R}. 

Similarly one can show that u > —ip in 11^(0) n < R}. Now we conclude by evaluating \(p\ in the 
region n^(0) n {|a;| < 1}. 

□ 

The following holds as a corollary of above lemma with 5^0: 
Corollary 2.10. Let tu satisfy the following in the viscosity sense: 
(a) -V+{D^uj) < 0, -V-{D^uj) >Q mT.R:= n^(0) n {|a;'| < R}. 
(h) duj/dv = {) onrjvnS^; 

(c) Lo — onVijC] Sfl; 

(d) \uj\ < R^-'. 
Then we have 

\uj\ < CR-" in n^(o) n {\x'\ < 1}. 

From Corollay |2.10 the following holds. 

Theorem 2.11 (general comparison). Suppose u,v : 11,^(0) — !■ M satisfies the following in the viscosity 
sense: 

(a) F{D^u,x) < < F{D^v,x) in U^{0); 

(b) u < V on T]j, \u\, \v\ < \x\'^^^ for large \x\; 

(c) % < fix) <lf, on Tn, 

where f{x) : M" — > M is continuous. Then u < v in 11,^(0). 
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2.4 Homogenization of the operator 

Here we state the homogenization resuh obtained in [8]. 

Theorem 2.12 (jH|). Let F satisfy (Fl) — {FA). Then there exists a unique operator F{D^u) : A^" ^• 
M" satisfying (Fl) — {F3) such that the following is true: 

Let he a hounded domain in M" whose houndary is and consists of two nonempty parts Td and 
Tn, and suppose that g,f : K" — > M are continuous. If u"^ solves 

' f) = 171 E; 

d^u^ = g(x) on Fjv; 

= f{x) on Yd 

Then uniformly converges to the solution u of 

' F{D^u) =0 in S, 



dyu = g(x) 
. u = f{x) 



on Fat; 
on 



Moreover there exists < a < 1 such that 



(9) 



Remarks 2.13. By Lemma 2. IF the theorem is valid in the setting of the unbounded domain f2 ~ ^u{0) 
with T]j — {x ■ V — —1}, Fat — {x ■ v — {)} and with uniformly hounded u^, though in this setting we 
would not have the quantitative estimate ([9|. 

Using above result we can apply compactness arguments to show the following. 



Theorem 2.14. Let K he a positive constant and let f 



he hounded and Holder continuous. 



For given v G 5" ^ , let : {~K < x ■ v < 0} 



he the unique hounded viscosity solution of 



(Pn). 



' F{D^UN,Nx)=0 m {-K<x-v<0} 
dvUN — f(x) on {x ■ V — 0}, 

u = 1 on {x ■ V = — -ftT} 



Then for any 5 > 0, there exists Nq only depending on K , the bound and Holder exponent of f such 
that 

\uN - u\{x) < 5 m {\x\ < K} forN>No, (10) 
where u is the unique hounded viscosity solution of 

F{D'^u)^0 m {-K < X ■ <0}; 

duU — f{x) on {x ■ u = 0}; 

n = 1 on {x ■ V = —K}. 
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Proof. 1. For fixed i', this is a consequence of Theorem 2.12 and Lemma 2.9 



2. Let us fix K and S. To show that iVo is independent of ly, suppose not. Then for some S > there 
exists a sequence of directions i^k and a sequence of uniformly bounded, uniformly Holder continuous 
functions fk € C'^(M") such that i^k t'o, fk{x) locally uniformly converges to f{x) and the the 



constants Nk = No{i/k) given by (10 1 goes to infinity as fc — )• oo. Let us denote u% and u'' by the 
corresponding solutions given in above theorem with v — Vk- Since {u'l,^} is uniformly Holder, along a 
subsequence it converges to a function uq. Using the stability properties of viscosity solutions as well 
as Theorem 1 2. 12[ one can check that uq solves 

F{D'^u) =0 in {-K < a; • i^o < 0}; 
< d^u — f{x) on {x ■ i/Q = 0}; 
u ^ 1 on {x ■ ly ^ ^K}- 



Now we choose A^o so that (10 1 holds for v = and 6/3, yielding a contradiction. 



□ 



3 Homogenization in half-spaces 

Now we are ready to study the averaging properties for solutions of (Pe), beginning with the strip 
domain: later we will approximate the general domain with these solutions. For given i> G 5", let us 
define 

H,(p) := {-1 < (x-p) -J. < 0}. (11) 

Consider u'^ solving 



{ i^(i?V,f) =0 in H,(p); 

lu = g{^-) on T^(v,p):^{{x-p)-v^Q}- 

u = \ on Ti{v,p) :— {{x — p) ■ v = —1}. 

The main theorem in this section is given below: 

Theorem 3.1. For given G S" and p G M", let solve (P^). Then the following holds: 

(i) For irrational directions v, there exists a unique constant /i(j^) G [mint;, max (/] independent on 
the choice of p such that locally uniformly converges to the linear profile 



u{x) — pl{v){{x — p) • + 1) + 1. 



(12) 



(ii) For rational directions v, ifV^ goes through the origin (that is if p — ^}, then the statement in (i) 
holds for V as well. 



(Hi) [Error estimate] Let u he an irrational direction. Then for solving (P^) and u as given in ( 12 I, 
we have the following estimate: for any < a < 1, there exists a constant C — Cq, > such that 



\u' - u| < Ce"/20 ^ Cu^{ef m n„(p), 



(13) 



where uj^ is as given in Lemma 2.2 (ii) 



10 



The proof of Theorem |3.1| consists of several lemmas below: the outline largely follows that of the 
corresponding result (Theorem 2.4) in |10| : for clarity we present the full proof below. 



Proof of Theorem [3TT1 



Due to the uniform Holder regularity of {ue} (Theorem 2.8 a)), along subsequences u^. — )■ w in 



Hi, (p) . Note that there may be different limits along different subsequences Sj . Below we will show that 
if v is an irrational direction then all subsequential limits of {u^} coincide. 

First we present a "flatness" result: we show that, sufficiently (e^/^^-) away from the Neumann 
boundary Tq = T(iy,p), is almost a constant on hypcrplancs parallel to Fq. 

Lemma 3.2. Let us fix v G 5" — MZ and p G M", and let us denote Fq = Fo(i^, p). Then for any 
xq € n^(p) with dist(a;o,Fo) > e^^^", and for any < a < 1, there exists a constant C = C{a,n) such 
that for any x G H(xq) := {{x — Xq) ■ v = 0} 



\Ue{x) - Ue(xo)| < £{e), 



(14) 



where £{e) := Ce"^^^ + Cu)i,{e)^ with lUi, : [0,oo) — t- [0, oo) as given in Lemma 2.2 (ii). 



Proof. Let xi G H{xo). By Lemma 2.2 (ii), there exists y G M" such that jsi — j/| < Me^/^°, y—xo G eZ" 
and 

dist{y,H{xo)) < eu!^{e). (15) 

Let us compare u^{x) with Ue{x) := Ug{x — Xq + y) in the domain E = Tl,y{p) nll,y{p + xq — y). Due 
to Theorem 2.8 ( 15 1 and the fact that y — Xq G eZ" it follows that 



\di,Us{x) - gi 



X- xo + y 



\dtyu^{x) - g{-)\ < U!i,{e)" on Fo - |(xo - y) • 

e 



(16) 



Moreover, due to the Holder continuity of up to the Dirichlet boundary Ti{i>,p) we have 

\ue{x) -l\< {euji,{e) f on Fi(iy,p) + |(a;o - y) ■ v\v 
Lastly, due to the fact that 1 < g < 2, we have |u6(a;)| < 2 in n^(p). 



Putting together above estimates, we can apply Lemma 2.9 to v{x) := u^{x) ~ Us{x) with 5 ~ uji,{e)°' 
and R» 5 io obtain that 



\u^{x) — u^{x)\ < Cuj„{e)'' in \x — xq\ < 1. 

Now we can conclude, by observing 

\ue{xo) ~ u^{xi)\ < \u^{xo) - u^{y)\ + \u^{y) - u^{xi)\ 

< Ciu,{er + \u,{y) - u,ixi)\ 



(17) 



where the second inequality is d ue t o (17 1 with x — xq and the third inequality is due to the fact that 
\xi -y\< M£i/i° and Theorem [2!6l 

□ 
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Next we define the average linear profile for u'^: pick a point yo E {x ■ = and for each e > 
let us define by fij such that 



(18) 



By Lemma 3.2 and by the comparison principle (Theorem 2.5 1, we obtain that 

U<u,<Um{-l<{x-p)-v< -e^/^"}, 
where U and U_ are linear functions which satisfies f7 = C/ = lon{(a; — p)-i^ = — 1} and 
U = u{xo) + E{e), U = u{xa) - £{s) on {{x - p) ■ v ^ -£-1/^°.} 



From (19 1 and the definition of Vj the following estimate holds: For x e Il^{p), 

\u,{x)-v,ix)\<£ie) + Ce'^^°. 
this implies that {u^} uniformly converges to linear profiles along subsequences. 



(19) 



(20) 



Lemma 3.3. [Theorem 3.1, (i)] If v is an irrational direction, then there exists a unique jin which only 
depends on v such that uniformly converges to a linear profile u(x) = fi^{{x — p) • + 1) + 1. 



Proof. 1. Let < 77 < e be sufficiently small, let pi,P2 G M", and let 



solve {Pj[) in Il„{pi) and 



u-q = Uri.p2 solve (P^) in Il^{q). We will show that \fi^ — fin\ can be made arbitrarily small if we choose e 
sufficiently small, independent of the choice of pi and p2- This concludes Theorem 3.1 (i). In the proof 
below we strongly use the fact that F(0, a;) = and thus that linear functions solve F{D'^u, x) = 0. 



2. Let us re-scale 



u^{£x) u^i-qx) 
e ?7 



and denote by Fi = ^To{i',pi) and F2 = ^To{i/,p2) as the corresponding Neumann boundaries of 



We 

and Wjf. We first translate Fi by r S Z" such that Fi and F2 are close. By (iii) of Lemma 2.2 there 
exist (71 G Fi such that 

P2 

— =qi+T + yi, 
V 

where \yi \ < rj and r G Z". Hence after translating a;^ by r, we may suppose that w^i^x) and Wn{x) are 
defined, respectively, on the extended strips 

rig := {a; : — < (a; — yi) ■ v <Gi\ and 17,, := {x : — < (a; — 2/2) • v < 0}, 
e r] 

where y2 = ^ and \yi - 2/2I < ?7- 

Without loss of generality, we may assume that y2 ■ ^ yi ■ v. Since \yi — 2/2I ^ and g G C^, we 
have 

\g{x)-g{x + yi^y2)\<v^. (21) 
From (21 ) and the C'^'" regularity of w^ix) (see Theorem 2.8), we conclude that there exists a constant 



Co > independent of e and 77 such that 

\d,w,{x)-g{x)\<Co{^f +7^^ 



on F„ 



{{x-y2)-y^Q}. 



(22) 
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Figure 2 



Let Ve be given by (18). Then by (20) 



(23) 



From (23) and the comparison principle, it follows that 



{fie-S{e)){{x-yi)-iy+-)<w,{x)-- <{fi,+£{s)){{x-yi)-,^+-) in ft,. (24) 

e e e 



(24) means that the slope of in the direction of ly (i.e. v ■ Dwe) is between /x^ ± £{e) on {x : 
{x — yi) ■ v = —^}. Now let us consider the linear profiles 

/i(a;) = ai{x - yi) ■ ly + hi and ^2(2^) = a2{x - yi) ■ v + 62, 

whose respective slopes are ai = /i^ +£{e) and 02 = /i^ —5(e). Here bi and 62 are chosen to match the 
Dirichlet boundary data of cj^. 



3. Now we define 

w{x) 



li{x) in {-l/v - V < - yi) ■ V < -l/e} 

We{x) + ci in { — 1/e < {x — yi) ■ v < 0} 



and 



w(a;) 



^2(2;) 



in { — — fy<(a; — zi)-i^<— 



We(x) + C2 in { — 1/e < (a; — j/i) • < 0} 
where ci and C2 are constants satisfying 

h = + ci, I2 = We + C2 on {{x - yi) ■ V = 

(See Figure 2.) 

Note that, due to (|24]), in {-i < (x - zi) • < 0} we have 

w{x) — Ta\a.{li(x),We{x) + ci) and u;(a;) = max(Z2(x), ^^(x) + C2), 
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and thus it follows that w and w satisfies, in the viscosity sense, 



F{D^w, -) < < F{D^w, -) 

e £ 



4. Let us define 



h,{x) = Co(ry« + rj^){{x ~ y2) ■ v + I/77), 



where Cq is given in (22). Then w 



= w + hi and w :— w~ hi respectively solves 
F{D'^w+,x)>0 in n^; 



di 



on 



and 



F{D'^w-,x)<Q in 
dw 



di 



on r^. 



Since 



~ Wri on {{x — p) ■ ly — —-}, from the comparison principle for (P^) (Theorem 



it follows that 

Hence we conclude 



2.11 



w < Wjj < in n.^. 



where /x^ is the slope of f^, which is defined as in (|18[) for m^. Now we can conclude. 



(25) 
(26) 

□ 



The proof of the following lemma is immediate from Lemma 3.3 and (26 1 



Lemma 3.4. [Error estimate: Theorem 3.1 (in)] For any irrational direction v there is a unique 
homogenized slope iJ.{v) S K and eg = ^oW) > such that for < e < £0 ihe following holds: for any 
p g M" and solving {P^) in 11^(0), 



\u,{x) - {pL{y){{x -p)-iy + l) + l)\< £{e) in n,(p). 



(27) 



Lemma 3.5. [Theorem 3.1 (ii)J 



Let f be a rational direction. If the Neumann boundary Fp p asses through p = 0, then there is a 
unique homogenized slope p.{i') for which the result of Lemma 3.4 holds with £{e) — Ce"/^. 



3.3 



Proof. The proof is parallel to that of Lemma 3.3 Let cj^ and be as given in the proof of Lemma 



Note that, since Sl^ and il,, have their Neumann boundaries passing through the origin, dwe/dv = 
g[xj = dwjjdv without translation of the x variable, and thus we do not need to use the properties 
of hyperplanes with an irrational normal (Lemma |2.2[ iii)) to estimate the error between the shifted 
Neumann boundary datas. 

□ 

As mentioned in [TU], if is a rational direction with p ^ 0, the values of g{-/s) on 957^ and Sfi,, 



may be very different under any translation, and thus the proof of Lemma 3.3 fails. In this case may 
converge to solutions of different Neumann boundary data depending on the subsequences. 
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4 Continuity of the homogenized slope 



In the previous section we have shown that for an irrational direction v e 5"^^ — MZ", there is a unique 
homogenized slope /i(i^) for any solutions of {P^) in Hv{p). In this section we investigate the continuity 
properties of ^{v) as well as the mode of convergence for u'^ as the normal direction v of the domain 
varies. For section 4 and 5 we assume the following additional condition on the homogenized operator 



F as given in Theorem 2.12 



F{M) : M' 



only depends on the eigenvalues of M. 



(28) 



As mentioned before, the condition (28) is equivalent to saying that F is rotation and reflection 
invariant. 



Theorem 4.1. Let ^i{v) : (5" 



be as given in Theorem 3., and suppose that F satisfies 



(28). Then fj, has a continuous extension ^{v) : S" 



More precisely for any v S''^ and 5 > 



there exists Eq = such that the following holds: 

If vi and V2 are irrational such that 



< 1^1 - \v2 -v\< £o, 



(29) 



then we have 

(a) \^i{vi) - n{v2)\ < 5; 

(b) For given p G M", the solutions u'^' of {P^') in 11^. (p) and the average slope /i(u^') given as in 
( 18 ) satisfies 



Remarks 4.2. In the proof we indeed show that, for any directions vi and satisfying (29), the range 



o/ {/u(u^')}£^i fluctuate only by S, if e is sufficiently small (see (60)). The fact that Vi's are irrational 
is only used to guarantee that there is only one suhsequential limit for ^(u^'), with i = 1,2. Note that 
our statement uses v as the reference direction and thus is a different type of estimate than those in 



Theorem 3. 1 ] 



Thanks to Lemma |3.4[ it is enough to consider the case p = Q. The main idea in the proof of 
Theorem 4.1 is to use u as the reference direction and approximate .g(-) by piecewise continuous 
functions, where each continuous parts are "projections" oi g on {x ■ u = 0} (see further description on 
the approximation below). For simplicity of the presentation, we will prove the theorem in M^: at the 
end of the proof will describe the modification required for M". 



4.1 Description of the perturbation of boundary data and a sketch of the 
proof 

First let us describe the main ideas in the proof. We begin by introducing several notations. For 
notational simplicity and clarity in the proof, we assume that v ~ e2'- we will explain in the paragraph 
below how to modify the notations and the proof for v ^ e2- Let us define 

r^o := nj.(0) = {a; e : -1 < ^2 := X • 62 < 0} 

and for z = 1 , 2 

n^.(O) = {a: e : -1 < X • < 0}. 
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Let us also define the family of functions 

gi{xi,X2) = gi{xi) = g{xi,S{i - 1)), where i = 1, ...,m := [\] + 1. (30) 

d 

(see Figure 3). Then gi is a 1-periodic function with respect to xi. 
oA remark for v ^ 62- 

In two dimensions, if is a rational direction different from 62, take the smallest e N such that 
K^v = mod Z^. Then we define gi{x) = 17(2;' + S{i — where x' = x — x ■ u, and gi is a _ftri,-periodic 
function. If v is an irrational direction, take the smallest K,^ G N such that IKj^vl < S mod Z^. Then 
gi as defined above is almost K^,- periodic up to the order of 5 with respect to x' . We point out that 
it does not make any difference in the proof divided in the following two subsections if we replace the 
periodicity of gi by the fact that gi's are periodic up to the order 6. 



o Proof by heuristics 

Since the domains f2i and Q2 point toward different directions and 1/2, we cannot directly compare 
their boundary data, even if dili and dfl2 cover most part of the unit cell in M"/Z". To overcome this 
difficulty we perform a multi-scale homogenization as follows. 

First we consider the functions gi (i — 1, .., m), whose profiles cover most values of g up to the order 
of 6^ , where P is the Holder exponent of g. Note that most values of g are taken on dfli and on dV,2 
since vi and i>2 are irrational directions. On the other hand, since i>i and 1^2 are very close to 62 which 
is a rational direction, the averaging behavior of a solution in fli (or 572) would appear only after e 
gets very small, as vi (or 1^2) approaches v = 62- 

Let N = [S/liyi — If — z/| = jz^i — 62] is chosen much smaller than S, then we can say that 
the Neumann data gi{-/e) is (almost) repeated N times on Tq = {x ■ vi = 0} with period e, up to the 
error 0{5^). (See Figure 4.) Similarly, on the next piece of the boundary, (72 ('Z^) (almost) repeated 
N times and then gzi'/e) is repeated N times: this pattern will repeat with gk (fc e N mod m). 

If N is sufficiently large, i.e., if \vi — v\ is sufficiently small compared to 5, the solution of {P^^) 
in f2i will exhibit averaging behavior, iVe-away from Fq. More precisely, on the hyperplane H located 
Nrj-Suway from Fq, Ue would be homogenized by the repeating profiles of gi (for some fixed i) with an 
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error of 0{5^). This is the first homogenization of near the boundary of fJi: we denote, by n{gi), 
the corresponding values of the homogenized slopes of on H in each A^e-segment. 

Now more than Ns away from dfli, we obtain the second homogenization of u^, whose slope is 
determined by n{gi), i — \, ..,m. In the proof this second homogenization is divided into two parts, 
in the middle region which is Ne to KNe-awsLy from Tq and then in the rest of the domain fli. The 
homogenization argument in the middle region is to ensure that the oscillation of the operator F in 
x-variable does not alter the behavior of the solution too much as v varies. Note that, due to ( [28| after 
rotation and reflection we may assume that the arrangement of ...fj,{gm) is the same for vi and 

U2. Therefore the second homogenization procedure applied to vi and 1^2 yields that — ^(^^2)! is 

small. 

Below we will present a rigorous proof for Theorem |4.H with above heuristics in mind. 



4.2 Estimates on localized barriers 

Let (5 > be given. We will first prove a series of lemmas which analyzes localized versions of in the 
three areas of the domain (near-strip, middle-strip, and inner-strip). Based on the estimates, in Section 

4.3 we will then put together the localized barriers to construct appropriate test functions (sub- and 
supersolutions) of (Pp). The goal is, to compare these test functions with u'^' to show the following: 



for given S > and for Vi and e satisfying the assumptions given in Theorem 4.1 there exists a constant 
/io such that 

Ia^«')-MoI <'5forz-l,2. 
From above inequality Theorem |4.1| follows. 

Let a 
Theorem 



and gi as given above. We continue with the notations. Given 5 > 0, let A^o be as given in 



2.14 



for Neumann boundary data / G with max |/| < max \g\. 
Let us choose irrational unit vectors 1^1,^2 G such that 

0<Vo<£o< mm[S^°JN^\ (31) 
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where Sq ^ {vi — and 770 = | — 62 1 . Let us also define 



and lastly suppose 



n ^ ^ s: 21/20 



(32) 



(33) 



The condition (33) is to ensure that 2mNe « S, i.e. we consider sufficiently small e so that to 
approach the correct averaged slope ^{vi). 

With above definition of e and N , consider the strip regions 

Ik = [{k - l)A^e, kNe] x M for fc e Z. 

Let k e denote k in modulo to, where m — [-] + 1. Note that, since \vi — 62] << 5, the function 

5 



gk{xi/e) defined in (30) is (almost) repeated TV times on 1^ H Fjv, where Fat :— {x ■ vi — 0}. This fact 



and the Holder continuity of g yield that 



-) - < CS^ on Fi n h for fc G Z. 

£ e £ 



(34) 



Similarly one can define Ik '■— [(fc — l)Af£, kMe] x M for fc e Z to approximate g on {x ■ V2 = 0}. 



4.2.1 Estimates for solutions near the Neumann boundary 



Let us first consider the averaging behavior of solutions near {Ne- close to) the Neumann boundary. 
For the strip domain := {— 27V£ < x ■ vi < 0}, let solve 



{N) 



( F{D^u,, x/e) = in {-2iV£ < x • J^i < 0} 
^(a^)=5(-) on {x-vi^{)} 

on {x-vi = -2N£} 



I We = 1 

Similarly we define iig in the strip domain 

% {-2Af£ < a; • :^2 < 0}. 



First we replace the Neumann data g with the family of functions gk introduced in ( 30 ) 
Let We and We be solutions in and H^, satisfying 

f F{D^We,x/e) =0 in {-27V£ < x ■ i^i < 0} 

—-^{x) = gk{ — ) on {x ■ i/i — 0} O Ik for fc e Z 

UVi £ 



. We = 1 



on {x- 1^1 = ~2Ne} 
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and 



' F{D^We,x/e) = in {-2Me <x-U2<0} 

-7r^ix)=9k{ — ) on {a; • z/2 = 0} n /fe for fc e Z 

, = 1 on {x ■ 1^2 = —2Me}. 

Next we localize the Neumann boundary data: for given k & Z, let Wg^i- and Wg^k solve 

( F{D^Ws,k,x/e) = in {-2Ns < x ■ i/i < 0} 
-^^{2:) =5ft(y) on {a; •1/1=0} 



and 



. We,k = 1 on {a; • z/i = — 2A''e} 

f F{D'^We,k, x/e) = in {-2M£ < a; • z/2 < 0} 

^ {x)=gk{ — ) on |x-:/2=0| 

oz/2 e 



on {x ■ U2 = -2Me} 



Wc will show that the profile of Ug and ilg does not change too much with above perturbation and 

localization of the Neumann data. 

Lemma 4.3. There exists a dimensional constant C > such that the following holds: For each k G Z", 
\u^-Ws\, \ws-w^^k\<C5'^Ne on {x ■ vi = -Ne] r\Ik 

and 

\ue-We\, |w£ - We.fel < C(5'^/^Me on {x ■ V2 = -M£}r\ik. 
Proof. We will only prove the lemma for Ue and Wg. Let us define 

Td ■■= {x- vx = -2Ne}, H := {x ■ ui = -Ns}, V n := {x ■ Vi = Q}. 



1. Since = on To, and 



on To, wc obtain that \ti^ — < 25^ Ne by adding linear profiles ±5^{x ■ ui + 2Ne) to u^, and to We 
and then apply the comparison principle to get respective orders. 

2. Next we compare and w^k- Let p G //j n Fq. Observe that 

\gk - gk+i \ < 5^ 

by the construction of gk and by the Holder continuity of g. it follows that 

\9k+i - 9k\ < 5^'^ for - 5-^'^ < I < 5-^/\ (35) 

Above inequality implies that 

\du-,Ws - dv^Ws^k\ < ^^^'^ in -B^-zs/^jveW ^{x-vx= 0}. 
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Lastly note that, since 1 < g < 2, we have 

\we - w,,k\ < 2A^e on dBaip) n 

to Wr 



We would like to apply the re-scaled (replacing Ne as the unit scale) version of Lemma 2.9 
and We^fc in Bii{p) n with R = S~^/^Ne. We then obtain 

\we - We,k\ < Cd^'^Ne on H n 4. 

□ 

Next we show that Wi^^k and We.fc are almost linear, and that they do not differ from each other very 
much. Let iiiw^^k) be given as the approximating slope of w^^k as in ( 18 1, i.e. choose a point pk € Hnik 
and define 4>e^k by 

(f)e,k = /i(We,fc)(2: • 1^1 + 2iVe) + 1, n{WeM) = ■ 

Similarly one can define (f>^^k and ^[we^k)- 



Lemma 4.4. Let a ~ a{n,K,\) he as given in Theorem 2.7 and let a\ = min[a,/3]. Then there exists 
a dimensional constant C > such that the following holds true: Respectively on {x ■ vi = —Ne} and 
{x ■ V2 = -Me}, 

\we,k - (t'sA < CNed"^ and \we,k - ~4>e,k\ < CMe6"K (36) 

Moreover 

\^^{w,^k) - Kwe,k)\ ^ CS'''^\ (37) 

Proof. 1. We will prove the first inequality on w^^k- Fix xa E H . Due to the fact that — ei| < 5, the 
following holds: for any y E H, there exists z such that 

Xq = z mod eZ'^ , \zi--yi\<e, d\st{z, H) < 6e. (38) 

Consider the function 

Ipe.kix) = We.k{x - (xq - z)), 

which still solves F(D'^ijj, |) = in lie + (^^o — z). We will compare w^.k and ■ipe,k in the domain 

I]:=n,n(n, + (a;o-z)). 
Without loss of generality, we may assume 

E = {-2Ne + Soe <x-i^i< 0}, 
where < Sq < \{xo — z) ■ vi \ < 5. Observe that comparison with linear profiles yields that 

\wE.k - ■0£,fe| < (^e on {x -vi^ -~'2Ne + Soe}. 



Moreover, since gk{x) is constant in a;2-variable, the last property of (38) and Theorem 2.8 yields that 

\d^^w^^k ~ d^ii^^^k] < C5°'^ on {x ■ vi ^ 0}. 
Let h — C5°'{x ■ vi + 2Ne) + 5e. Then by applying the comparison principle in E we have 



,k < tpe,k + h and ^e,k < We,k + ft. in E. 
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Since <h< 2Ne5°' on H := {x ■ vi = -Ne}, we get 

\w£^k{xo) - We,k{z)\ = \we,k{xQ) - V'e.fe (a;o) | < 27V£(5"i on H. 
Now we conclude that 

|w£,fc(xo) - We,fc(2/)| < \We,k{xo) - 'We,k{z)\ + \We,k{z) - We,k{v)\ 

< 2Ne5''' + \weMz) - WeMy)\ 

< 2NeS°'' + C\z - y\ 

where the last inequality follows since NS" ^ S^'^" /eq > 1. 



(39) 



3. To prove the statement for Ws, one can argue as above, by replacing Ik by Ik and by Ilg, and 
H with H := {x-i^2 = -Me}. 



4. We proceed to prove (37). Recall that AI > N. First we compare w^^k with p^, solving 
' F{D^p£, f ) = in {-2Ne < a; • < 0}, 
d„2Pe = 9k{j) on {x- 1/2 = 0}, 



We will show that 



Pe = 1 on {a; • = —2Ne} 



,k\ < CS"^^'^Ne in {-2Ne < x ■ V2 < 0}. 



(40) 



If (40) holds, then we can conclude by barrier arguments as given in the proof of Lemma 3.3 to show 
that and Co^ are close. 

Due to (36), it is enough to argue in the iVe-neighborhood of the origin. Note that the Neumann 



boundary of p^ and ^ only differ by R5e in i?A^e-neighborhood of the origin. Therefore Theorem 2. 
yields that 



\du^Pe - 9k{-)\ < {RST' on {x-i/,=0}n {\x\ < RNe}. 
e 



(41) 



that 



Moreover, since 1 < g < 2, we have jp^ — We^fej < 2Ne on {x ■ vi = —2Ne}. Therefore (41) yields 
\pe-w,,k\<C{RSr'Ne on{x-i^i = -Ne}nlk+j,\j\<R. (42) 



Now due to (41) and (42), Lemma 2.9 yields that 

\Pe~w,^k\<C{{R5r' +R-^)Ne. 
Let us take R = 6^^^^ to conclude. 



□ 



The following Corollary is immediate from Lemma 4.3 and Lemma 4.4 as well as Theorem 2. 



Corollary 4.5. Denote pk f^{we,k), then for ai given as in Lemma 

, d 



u,-pk\ <'5"i onHnlk. 
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4.2.2 Estimates for solutions in the middle region 

From the previous section we have seen that the solutions and Ue average in the unit of Ns, when the 
solution is TVe-away from the Neumann boundary, and moreover that they average to a similar value. 
If the operator F was homogeneous we could use this result and perform a second homogenization, and 
use (28) and (37 1 to conclude. However we have to be careful with the inhomogeneities of F, such that 
different normal directions of the hyperplanes pointing toward i^i and 1^2 does not make a difference in 
the way the respective solutions average out. This is what we are going to analyze in this subsection. 



Let K 



J, and let us consider the domain 



E = {a; e : -2KmNe < x ■ i^i < ~Ne}. 



Note that by (33) 2KmNe < 1. Let fik as given in Corollary |4.5[ and let solve 

F{D^v,,-/e) 



dVe 



in S 

on _ff n Ik ioT k G Z 
on {x ■ vi — —2KmNe}. 



Next let Ve solve 



F{D^v,) = 



1 



on Hnik for A: e Z 



on {x ■ v\ = —2KinNe}. 



We will show that and are close to each other. 
Lemma 4.6. There exists a — a{n,A,X) such that 

\Ve{x) - Veix)\ < S'^NS 



in E. 



(43) 



Proof. Let us apply Theorem 2.14 with f{x) — jik on H (1 Ik (or its interpolation so that f{Nex) is 
Holder continuous) to the re-scaled function 



Due to (jSlj), ^ > iVo, and thus Theorem |2.14| yields ( |44[ ). 

\ve{x) -Ve\< SNe in Ei. 



(44) 

□ 



Next we are going to show that is close to a linear profile away from H. To see this, recall that 
the slopes fii,...,firn are repeated if-times on the Neumann boundary of E. Hence the homogenization 
arguments for the homogenized operator F should apply in our setting, if K is chosen sufhciently large. 
This is what we will show below in detail. 



Let fJ.{ve) and ^(we) be respectively the average slope of and respectively in E as in (18). 
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Lemma 4.7. Let a given as in Lemma 4-4 Then there exists a dimensional constant C > such that 



d 

I— - ^i{ve)\ < CS" on L^ {~2KmNe < x ■ vi < -KmNe}. 
ov\ 

Proof. Choose a point xo G L. Since has constant slope fik on H D Ik with \H D Ik] = Ne, and since 
the slopes fj.i,...,fim are repeated on H, for any y such that (y — xq) ■ vi ~ there is a point z e L such 
that 

\y - z\< mNe, xq - z e H, (45) 

and 

d _ d _ 

—v^{x) = -w—v^{x - {xo - z)) on H. (46) 



Due to (46) and Theorem 2.5 we have v^{x) = v^{x — {xq ~ z)). Moreover 

\veixo) - Ve{y)\ = \ve{^) " vM\ < CmNs (47) 
wher e the i nequ ality follows from the interior Lipschitz regularity of xmNe {KfnNex) (Theo- 



rem 



2.7|). With (|47|), we can conclude due to the C^'" regularity of ^ 



KmNe 



v^{KmNex). 



□ 



Similarly, for the normal direction 1^2 one can construct and accordingly in 

E = {-2KmMe < x ■ V2 < -Me}. 
Note that Corollary |4 . 5 1 implies 



(48) 



Note that, due to (28) after rotation and reflection we may assume that the arrangement of /i((?i) is the 

(49) 



same for vi and V2. Therefore from (48) and (28) it follows that 



IM^^s)-/i(^s)| <^"- 



Due to Lemma 4.6 there exists a = a(n, A, A, /3) such that 

HVe)-^^iiie)\<S''. 

4.3 The proof of the main theorem 



(50) 



Recall that u^' for i = 1,2 solve (Pp) in ll^{0). With the estimates obtained in the previous section, 
we are now ready to construct our barrier for u^^ and . As mentioned before we will construct 
barriers in three separate regions and patch them up. Let us construct it for the normal direction vi: 
let ao = min(a, /3/4) where a is as given in (50). 

In far-away region: let us define 
(F) f,ix) ■.= A{x- 1^1 + 1) + 1, 

where 

A = fi{v') + 106°'°. 
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In the middle strip: next consider p^: the unique (bounded) viscosity solution of 
' Fp^^e, f ) = in {-KmNe<x-vi<-2Ne}; 

(M) J |£^^A2 on i7 = {x • i/i = -2iV£}; 

— fi, on {x ■ vi ~ —KmNs} 

where A2{x) is obtained by approximating pk such that 

A2(x) e C^(M") with its norm less than 5'^{Ne)''^; 



and 



A2 is periodic with period mNe on H, and /i/c + 26^ < A2{x) < pk + 5(5^. 



Such A2 satisfying (51)-(52) exists for ao < /3/4 since, due to (35l, for any given fc,^ e Z, 

\pk+i ~ Pk\ < |^|<5''- 

In the near-boundary region: last, let (p^ be the unique (bounded) viscosity solution of 
' F{D^(t>e,f)=0 in {-2Ne < X ■ vi < -2Ne}; 

(N) J |^=ff(f) on {x-i.,=0} 

4>e = Pe on {x ■ vi — —2Ne}] 

Our goal is to show that the function 

/e in {-1 <x-vi< -KmNe} 

Ue := < Pe in {-KmNe < x ■ vi < -2Ne} 

in {-2Ne < x ■ vi < 0} 

is a supersolution of (P^) in H^^. 



(51) 
(52) 

(53) 



Observe that, due to Lemma |4.6| and Lemma |4. 7 

dp, 



dvi 



< K on {x ■ vi — ~Km,Ns}. 



(54) 



This means the patch-up of and together is a supersolution, i.e., 

F{D^U^, ^)>0 in {-1 < a; • J^i < -KrnNe}. 
It remains to see if the same holds for the patch- up of and 

Note that p^ is a constant on its Dirichlet boundary {x ■ i^i = —KmNe}. Due to this fact and the 
small oscillation of A2 given in (51), as well as Theorem 2.8 and Theorem 2.14 applied to re-scaled 
versions of p^ yields that -^p^{Nex) has its C^'^ norm of size 6^^^ in the scale of Ne . More precisely, 
for any y E H := {x ■ i^i = —2Ne} we have 



\pe{x) - pe{y) -IJO-ix- Xo)\ < (5"°iV£ On H D B25-.'o Neiv) , 



(55) 
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where 770 = Dps{xo) - Dp^{xo) ■ vi. 

Let us fix a point G H. Next we will invoke Lemma 2.9 for ■^^/^{Nex) with A 
R = 26-"° to obtain 

\(f>e -m-{x- xo) - 0e| < S''°Ne in {-2iV£ < x ■ i^i < 0} D Bg-.o^M, 
where 4>e solves (N) with the fixed boundary data pe replaced by the constant pi,{xo). 
Suppose xq £ H O Ik- Due to Corollary |4.5| we have 



\p{(t>e) - /ifcl < ^""Ne on H -.^{x-vi^ -Ne}. 



(56 1 as well as (57) yields that 



6°"° and 



(56) 



(57) 



\(f)e{x) - {rio ■ {x - xo) + pkix - a;o) ■ ^^i - (t){xo)\ < 2S°'°Ne on {-2N£ < x ■ vi < ~Ne}. 



It follows from (551,(58) and the C^'" regularity of and 0^ in TVe-scale that 

d , d 

-^Ve < ^Pe on H, 

and thus the patch-up of p^ and (j)^ is a supersolution, i.e. 

F{DU'', ^) > in {-KmNe < x ■ vi < -Ne}. 



(58) 



Summarizing, we have shown that is a supersolution of {P^'^ ) in (0) . Hence by comparison 
principle we obtain that < in In particular. 



<i [x) < Ue{x) < {^i{v^) + 10(5"°)(x • - 1) + 1. 
Similarly one can construct a subsolution Ve, to show that 

\u'','{x)- p{v^){x-iyi - 1) - 1| < 10(5"°. 

Above equation yields that 



(59) 



Similarly, 



Moreover by ( 50 1 we have 



|m«^)-M^^^)|<105"°. 
|m«^)-m(5'')| < 10S"«. 
Hvn-^^m\<s?. 



This proves Theorem 4.1 (a). 

Lastly observe that {P^'^) given in section 4.2.2, has its only microscopic term in the Neumann 



boundary data, which are constants in each segments Ik of H of size 

constants repeating itself after j pieces of Ik- that is the Neumann boundary data is periodic with period 
£^/^°. Since solves (Pp), fJ-{ve) converges uniformly to a limit /io as e — with the convergence 



rate only depending on S and e, as long as ( |33[ ) is satisfied. Therefore, if necessary approximating the 
Neumann data with functions of the form /( 1/20 ) where / is periodic and Lipschitz, we conclude that 

\^i{u'',') - Pol < ioT I ^ 1,2, (60) 

if £ is sufficiently small. Since (5 > is arbitrarily chosen and ao only depends on f3,n,X and A, this 



proves the second claim of Theorem 4.1 (b) 



□ 
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Remarks 4.8 (For dimensions higher than two). In higher dimensions n > 2 and for v = e„, one can 

define 

gi{xi,...,Xn-l,Xn) = gi{xi,...,Xn-l) ^ gixi, Xn-1, S{i - 1)) 

for z = 0, 1, rn = [S^^] . Let us also define 

4i,fc2,..../c„_i = [{ki - l)Ne,kiNe] x ... x [(A:„_i - l)Ne,k„^iNs] x M. 

Then parallel arguments in the previous two subsections apply to yield the corresponding result to The- 
orem \jn\ 



5 In general domain 



In this section we will use Theorem 4.1 as well as stability properties of viscosity solutions to prove our 
main result. As given in the introduction, let f2 be a bounded domain with boundary containing a 
unit ball K — _Bi(0). Let us denote v = the outward normal vector of f2 at x € dCl. Suppose that 
90 does not have any flat boundary parts in the following sense: For any xq € and sufficiently small 
a, there exists > and r{a) > such that 



Wx — i^xol > cr if a; e dfl n {r{(j) < \x — xq\ < ro}. 



(61) 



Note that, for example, any strictly convex domain satisfies (a)-(b). We now state the main 
theorem. 



Theorem 5.1. Let fl and K as given above, and let fx : 5" — > 

Consider solving 

' F{D^u^, 2^) = in n-K; 



1,2] be as given in Theorem 4-1 



u = 1 



on dfl; 
on K, 



and let F be the homogenized operator given in Theorem 
converges uniformly to the unique viscosity solution u of 



2.12 



with the assumption ( 28 1 . Then u' 



' F{D'^u) =0 m n-K; 
p.{v) on dfl] 
u = 1 on K. 



du 
du 



Remarks 5.2. Note that the uniqueness of u follows from the continuity o//i(i^) and Theorem 2.5 



Next we locally approximate with the solutions associated with strip domains discussed in section 
4, based on the regularity properties of and dQ. 

For given domain il , p e dfl and < fc < 1, let us define — Vp and 
Sfc := n {x : -e'' < {x ~ p) ■ i^o} D B^5k/s{p). 
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Let Wp solve 



f F{D^w„^) = in Efc 

on dfl n aSfc 
on 9Efe — c?r2. 



dv 



Wf = 1 



Next let 

and let solve 



Sfc := {-e'' <{x-p)-vo< 0} n S^s^/s (p), 



1 



= in Efe; 
) on {(a; -p) • i^o 0}; 

on {—e'' = V ■ [x — p)}. 



Since dVl is C^, we may assume that the hyperplane {{x — p) ■ vq — 0} is contained in the e^^/^- 
neighborhood of dVl in B^5k/s{p). 

Lemma 5.3. There exists < fc < 1 and a > 0: independent of e,p and vq such that 



vJ <e 



k-\-a 



in Efc n Bp2k/3{p). 



Proof. After a translation, we may assume that p = 0. First note that We and will oscillate at most 
of order in their respective domains and S^: This can be checked by comparison with linear 
profiles. Let us consider the re-scaled functions 

'w{x) — Wp{ex)/e and v{x) = v^{ex)/e. 



Then Theorem 2.8 the Holder continuity of g, and the fact that w and v oscillates up to Ce'' ^ yields 



in their respective domains -E^ and -Efe, where < a < 1 is as given in Theorem 
we have 

\d^,w, - d^,Ve\ < 0(e'=-i+(^-i)") on H := {x ■ ivo = -e^''/'^} n 5,5. 
Let us choose k sufficiently close to 1 so that 

,5k 



Consequently 
(62) 



Let us define 



k-l + {— - l)a> a/6. 



A 



h{x) := e-''/^\x - p\^ + n-£-'=/4(^2fe „ - p) . ^.^p), 

A 



-V^{D^h)>0, /i>e'= on EfenaSgfe/s, h > on {{x - p) ■ i^q = -e''} 



so that 



and ^ 0(e'"') << e"/^ on H, where H is given in ( [62| ). Now by (j62|, the comparison principle 
(Theorem 2.5) applies to oj, — and h in the domain E^ n {a; • z^o < —e^^*} to yield 
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\weix) - Veix)\ < Ce"^^ {{x - p) ■ i^a + e'') + h{x) in Sfc. (63) 



By evaluating the upper bound obtained in (631 in the region n B^2k/a, we conclude the lemma for 
a = min{Q;/6, fc/12}. 

□ 

We are now ready to show the main proposition. Let us define the semi-continuous limits 
limsup*M^(a;) ;= liminf(sup{u'^(j/) : y G H and |a: — ?/| < e}) 

and 

lim inf*u^(a;) := limsup(inf{it^(?/) : y G and |x — y| < e}). 

It is straightforward from the definition to check that limsup*u'^ is upper semicontinuous and 
liminf^u"^ is lower semicontinuous. 

Our main theorem is a consequence of the following proposition. 
Proposition 5.4. The following holds true: 
(a) u := limsup* is the viscosity suhsolution of (P); 
(h) u:= liminf* is the viscosity supersolution of (P). 

Proof. We will only show (a): parallel arguments apply to (b). The proof follows the perturbed test 
function method introduced by Evans |15| . and consists of several round of perturbations of the test 
functions (j). 

1. It follows from previously known results (see e.g., IS],[3) that 

F{D'^u) < in n. 

Moreover from comparison with linear functions it is straightforward to show that u < 1 on K. Therefore 
if u fails to be a subsolution of (P), then there exists a smooth function (p which touches u from above 
at a boundary point xq G dft and satisfies, for some S > 0, 

FiD^<j)){xo)>26 and |^(xo) > Ai(i^.J + 3,5 (64) 

Due to the regularity of and the continuity of /i(i^), there exists r > such that 

F(D'^4>)>S in r^nSJxo) and ^>fl(iy)+26 mdnnBrixo). 

ov 

Let denote the outer normal of dO. at Xq. By adding 5"^ {x — Xq) ■ + 5{x — xqY — M{{x — Xo)-vo)^ 
to 0(a;) and restricting the domain to Br{xo) with r < 6, we may assume that 4> > u except dX x — xq 
in f2 n Br{xo) and 

F{D^4>,^) >Omnr\Br{xo) and ^ > jlii^) + 26 in dn n Brixo). (65) 

2. Let us choose ao sufficiently small such that 
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|21/20 



> £ 



1-fe 



for X £ {Br{xQ) - Bern, {xq)) n dVl, 



(66) 



where < /c < 1 is the constant given in Lemma 5.3 This is possible due to our assumption (64), i.e., 
since d^l does not have any flat boundary parts. We now perturb and consider 0^ such that 

(a) 0^ uniformly converges to as e — > in fj n Br{xo); 

(b) F{D^',l)>OmBr{xo); 

(c) d^cf)" > - Ce"o on dQ n Br{xa), and d^cf)" > 2 on ar2 n 5^=0 (xq). 

Note that by (b)-(c) 0^ — can have its minimum at y in Br{xo) only if y is more than e^^-away 
from xq. This step is necessary to exclude the possibility that v^q is rational and the irrational normal 
direction near the minimum of 0^ — is too close to i^xo compared to the size of e to observe the 
averaging behavior. 

For convenience let xq = 0, v^a = e„, and set 7 := e"''. We will construct 0'^ as + (/^ + u^, where 
and will be constructed below. Let us first explicitly g = g-y, which satisfies g £ C^(i7), — >■ 
as 7 — > in il, and 

d^,g > -C-f'^ on dVL, d^g > 2 in B-^{xa) n dn. (67) 
Let us define S-^ {|a;'| < 27} x {—7 < x„ < 7}, and let us define g hy g = in fl — E^, and 



g{x',Xn) 



1 



where 0(r) : M [0,1] is a function which has support [—2,2], 0(r) = 1 for r G ( — 1,1) and 
0'(O) = 0. 

Note that, since dil is C^, dil n C {\xn\ < C-f^}: in particular, for small 7, 17 n {xn > 7} is 
empty in f2 n {|a;'| < 7}. It follows that g is in C^{Cl). Clearly g goes to zero in $7 as 7 —> 0, and so it 
remains to check ( 67 1 . Straightforward computations yield 

9.5(^)> ^(x„+7)M^)(l-C7)-C7(^||0'||oo(^(a;„+7)') > 2 on 5f7 n {|x| < 7} 

> -C7 in dfl, 



from which (67) follows. 

Next let V = Vj he the viscosity solution of 

-V+{D\,) = -V+iD-'g^) in fin B^ixo); 

d^v = on dnr]Br{xo); 

V — gj — on dBr{xo) n ft. 

To ensure that the remaining term v in the construction of cf)^ would vanish as e 
the LP norm of D^g. Observe that, from the definition of g, it follows that 

8 



0, we estimate 



Hence 



\g,j\ < - for i,j = 1, 
7 



(68) 
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and in particular if p = n, 



(69) 



Note that, due to (691, II/7IIL" is uniformly bounded. Therefore [H] yields that {v^}^^^ are uni- 
formly Holder continuous in Br{xo) n (Proposition 2.1 and Proposition 4.2 in [H]). Moreover, due 
to (68), 11/7 II LP vanishes as 7 — >■ for p < n. Due to this fact and the equi-continuity of v-y one can 
make use of the stabiHty of L^-viscosity solutions tor n > p > n — Eq for Eq = (n. A, A) (see [7]) and 
employ a compactness argument as in Lemma 2.3 of 20j to show that v-y uniformly converges to zero 
in f2 n -6^(2^0) as 7 — > 0. 

Due to above properties of g and v, (f)^ + ; 
note that 



Va with (J = e"'" satisfies (a)-(c): to check (b), 



-) - V-{D^g + D\) > F{D^(j), -) - V+iD^g) - V+{D\) > 0. 



4. Due to the definition of u and (a), for r > as given above we have > u'^ on dBr{xo) if e is 
sufficiently small, along a subsequence of e. Moreover, 0^ — has a local minimum at G Br{xo) D Cl 
with — >■ xo as e — >■ 0: we may add a constant to (j)"^ to assume that 4>'^{xe) = u'^^x^)- Let us denote 
fje — ■ Note that Xg lies outside of Bi^ao (xq) due to the construction of and due to the comparison 
principle applied to (j) a-nd u, G dfl D Bj.{xq). 

5. This step is to extract only the normal component of cj)^ at x = x^ and make a new function (p 
with it, so that we do not have to worry about the tangential derivative of 0^ at x^- We also adjust 
accordingly. Let us decompose ip'^ into 4>'^ ~ 4>i + 4>2 where 

Mx) ^{X- Xe) ■ [Dcj,' - Ve{% ■ Ddp')){x,). 
rji; ■ D(j)i{Xe) =Q, D(f)2{Xe) ^ VeiVe ■ D(j>2)iXe), a.Ild (f>2{Xe) ^ <p^ {Xe) ■ 



Then 



(70) 



Observe that, due to (70 1 and the fact that 0i is a linear function, 02 still satisfies (64 1 instead of 0. 
Furthermore, since 02 is smooth, we may choose e sufficiently small to replace 02 by a linear profile 

ip{x) := <j)2{Xe) + D(j)2{Xe) ■ [x - X^) 

such that 

02(2:) < <p{x) + Ce'^l^ in B,5./8(xe). 
where C depends on the C^-norm of 02- 

Next we will tilt accordingly so that we can use (/? instead of (\f as the test function. Let us define 

u,(x) = Ue{x) - 01 (a;) - Ce'^l^. 



Then satisfies 



Note that 



F{D'^u„x/e) = 

Ue{Xe) = (p{Xe) 



in D, 

in O n -B^sfc/s (xe) 
on dil n B^5k/s{xs). 



(71) 
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6. Next we will approximate il with a strip domain with the help of Lemma |5.3[ For e and as 
given above, let us pick e dfl such that 

(a) \pe - < e, 



(b) At dfl is normal to v^: an irrational vector. 
Let Ve solve the problem 

' F{D^Vs,x/e)^0 in ■= {-e'" < {x - Pe) ■ t^e < 0} 

lv^ = (p on {(a; -Pj) • = 



Due to Lemma 5.3 we have \lu^ — < e'=+" in 57 n Ee n Bg2k/3{x^), where solves 
' F{D^uje,^)^0 in nn{-e'' <{x-Xe)-iye}r\B^5k/s{xe), 

llf =5(f) + Ce''/^ on dnni^e" <{x-p,)-j^,}nB^..,s{x,), 

\^ OJe = f on {(x -Pe) • i^e = -e''} U 9B^5fc/8(a;e). 

Now comparison principle applied to and cj^ in {— e*^ < (a: — x^) ■ v} DilD B^sk/s{xe) yields that 

Ue <We <V, + in fl (1^,^(1 B^2k/3{Xe). 



Hence we have 



Combining (71) and (72), 



ip{x,)<v,{x,) + e''+'' + Ce"'/\ 



(72) 



(73) 



7. Lastly we will use Theorcm |4.1| to derive a contradiction to ( 73 ). Let us choose G eZ",|x£— j/gj < 
e and define 

Veix) := e'''ve{s''{x + y,)) 

so that Vg solves (P^i^^) in nj,e(ze) with jzej < e^^*^. Due to (|66| and Theorem 
sufficiently small so that 



4.1 



we can choose e is 



Then due to (65), the comparison principle applies to (p{x + y^) — -^{{x ~ z^) • + 1) and to yield 



that 



After scaling back to v^, we get 



,(0)<e-V(ys)-2- 



31 



Since \ve{xe) — Ve{ye)\ < Ce" for any a due to Theorem 2.8 we choose a > fc to obtam 



if e is sufficiently smah; this yields a contradiction to ( 73 ) , and we can conclude 



□ 



Proof of Theorem 15.1 



Due to Proposition 5.4 and the comparison principle (Theorem 2.5), we have u < u. But due to 
definition u > u. Therefore we conclude that u = u: this is equivalent to the local uniform convergence 
of the sequence to u ^ u ~ u, which is the unique viscosity solution of (P) . 

□ 
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